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simply pressure regulators, which produce aerodynamic noise when in 


Abstract — Piped natural gas is commonly distributed by regional gas 
distribution companies, which have pipelines that are responsible for 
distributing this fluid to final consumers. These meshes operate with 
maximum pressures in the order of 40 bar (4.0 MPa), however, there are 
other pressure levels involved, either for operational or safety reasons. The 
reduction of these working pressures is carried out in devices commonly 
called pressure reducing stations (EPRs), which have interconnected 
valves and accessories that enable the referred pressure reductions as 
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60534-8-4, Matlab, distribution networks. operation due to their physical configuration, local flow conditions and 


their own functionality. This paper presents a proposal for the 
determination of aerodynamic noise present in pressure regulating valves 
of piped natural gas distribution networks, based on the resulting 
mathematical formulation, predicted and recommended in IEC 60534-8-4 
“Prediction of noise generated by hydrodynamic flow” code. From these 
assumptions a computational code was developed in Matlab environment 
that allows the operations engineer to simulate and evaluate the noise level 


present in a certain operating condition, favoring his decision making. 








I. INTRODUCTION as required. These stations thus include pressure control 


Piped natural gas is commonly distributed through valves of various types (axial and globe, for example), 


regional gas distribution companies, which own pipelines depending on the flow rates involved and the design 


that are responsible for localized transportation and configuration considered by each distribution company. 


distribution of this fluid to end consumers, whether 
industrial, commercial, or residential, among others. This 
distribution meshes operate at lower pressures than 
transmission pipelines, reaching maximum pressures in the 
order of 40 bar (4.0 MPa), although there are usually still 
other pressure levels involved along the meshes, either for 
reasons of concern operational or safety. 


The reduction of these working pressures is performed 
in devices commonly referred to as pressure reducing 
stations (EPRs), which feature a set of interconnected 
valves and accessories that enable said pressure reductions 
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Pressure control valves or simply pressure regulators 
produce aerodynamic noise when in operation due to their 
physical configuration, local flow conditions and their own 
functionality. 


This noise has been, in many situations, a focal point of 
attention by regulatory and supervisory bodies, triggered in 
part by local legislation and current occupational health 
standards. Between 70 dBA and 85 dBA, the noise source 
becomes quite pronounced. However, long-term exposure 
to noise at this level is unlikely to damage human hearing. 
At 90 dBA and above, the noise source reaches a level 
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where sustained exposure can begin to damage human 
hearing. For these reasons, guidelines from the 
Occupational and Safety Hazard Administration (OSHA) 
mandate an 90 dBA limit for exposure of up to eight hours, 
as reproduced in “Table 01”. 


Table.1: OSHA Exposure time limits for various noise 
levels (Liptak, 2006) 


Duration of exposure (h) | Sound level (dBA) 
8 90 
4 95 
2 100 
1 105 
0.5 110 
0.25 or less 115 (max) 





The present work presents, therefore, a revision 
corresponding to the mathematical formulation of the 
noise generated during the operation of pressure control 
valves (according to IEC 60534-8-4 “Prediction of noise 
generated by hydrodynamic flow”), as well as the 
development of a mathematical algorithm, developed in 
Matlab environment, which allows a quick assessment of 
this problem by the operations engineer responsible for a 
certain piped natural gas distribution network. 


IL. MAJOR SOURCES OF NOISE 
2.1 — Mechanical vibration 


Mechanical vibration is caused by the response of 
internal components within a valve to turbulent flow 
through the valve. The noise generated by this type of 
vibration presents a worrying feature, as resonance 
problems may arise in the valve, a fact that would certainly 
cause fatigue failures in the internal parts of the valve. In 
practice, however, noise from mechanical vibration is 
infrequent in control valves. It is noteworthy that there is 
currently no reliable method for predicting noise generated 
by mechanical vibration in control valves. 


Mechanical vibration is generally below 100 dBA. 
2.2 — Hydrodynamic noise 


The Liquid flow noise, cavitation noise, and flashing 
noise can be generated by the flow of a liquid through a 
valve and piping system. Of the three noise sources, 
cavitation is the most serious because noise produced in 
this manner can be a sign that damage is occurring at some 
point in the valve or piping. 


Hydrodynamic noise caused by liquid turbulence, 
cavitation, or flashing is generally below 110 dBA. 
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2.3 — Aerodynamic noise 


Aerodynamic noise is a direct result of the conversion 
of the mechanical energy of the flow into acoustic energy 
as the fluid passes through the valve restriction. The 
proportionality of conversion is called acoustical 
efficiency and is related to valve pressure ratio and design. 
The aerodynamic noise is the most important component 
of the acoustic problem of a control valve, since it is 
generated by the pressure waves produced by the fluid 
turbulence or by other fluodynamic phenomena connected 
with supersonic waves. Aerodynamic noise can reach 150 
dBA. 


Aerodynamic noise prediction is based on the 
equations and nomenclature of the international standard 
for control valve noise prediction, IEC 60534-8-4. 


IHI. AERODYNAMIC NOISE PREDICTION 


The conversion of static pressure to kinetic energy at 
the vena contracta in the control valve creates high- 
velocity jets which can be subsonic, sonic or even 
supersonic. Turbulence and sonic shock waves create a 
noise spectrum with a characteristic peak frequency. 


The standard IEC 60534-8-3 - consists of a mix of 
thermodynamic and aerodynamic theory and empirical 
information. The design of the method allows a noise 
prediction for a valve based on the measurable geometry 


of the valve and the service conditions applied to the valve 
(Emerson, 2017). 


The method describes two different noise sources that 
can contribute to the overall noise generated by the valve: 
trim noise and valve outlet noise. The trim noise is 
dependent on the type of trim and its geometric features. 
The valve outlet noise is dependent on the valve outlet 
area, valve outlet Mach number and any expander 
downstream of the valve. 


The flow regime for a particular valve is determined 
from inlet pressure, downstream pressure, fluid physical 
data, and valve pressure recovery factor. Five flow 
regimes are defined as (ANSI/ISA, 1991): 


Regime I - Subsonic 


Regime II - Sonic with turbulent flow mixing 
(recompression) 


Regime III - No recompression but with flow shear 
mechanism 


Regime IV - Shock cell turbulent flow interaction 


Regime V - Constant acoustical efficiency (maximum 
noise) 
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The control valve aerodynamic noise prediction 
flowchart illustrated in “Fig. 1” illustrates the steps to be 
considered by the designer. 


Mass Throtiling jet 
flow diameter 


Jet velocity Based on vena 


contracta 
velocity 


calculation 


Stream 
power 


Pipe transmission 





Acoustical Sound power 


efficiency 





Peak 
freauency 


Internal sound 





loss pressure 


Fig. 1: Aerodynamic noise prediction flowchart on a 
control valve (Dresser, 2002 — adapted) 


} 


3.1 — Throttling jet diameter 


D; = 4.6(107*)F, CF (1) 
where Cy is the valve capacity and Dj is the equivalent jet 
diameter at the trim output, which is a typical constructive 
data for each trim type. It is directly proportional to the 
trim form factor Fa, whose typical values are easily found 
in the reference bibliographies (for a fully open circular 
flow orifice: Fa = 1). The most important of these 
parameters is the recovery coefficient F., which, at 
subsonic flow conditions, represents the energy fraction 
wasted inside the valve. 


3.2 Flow regime 


Regime I — Subsonic — this deals with all flow 
conditions in which the pressure drop is less than or just 
equal to the pressure drop that would produce sonic flow at 
the vena contracta. The flow is therefore subsonic and 
there is isentropic recompression. 


Py > Pe = Poc (2) 


where pzc is the downstream pressure, corresponding to 
the threshold of criticity. 


Poo = Py — F ip — Prec) (3) 
and 


a a 
Pree — PL r (4) 


where, pycc is absolute vena contracta pressure at critical 
flow conditions and y is the ratio of specific heats. 


Under this condition part of the mechanical energy 
existing in vena contracta is recovered as pressure energy 
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downstream the vena contracta. The remaining energy is 
wasted by turbulence, thus changing into heat and noise. 


Regime II — Sonic with turbulent flow mixing 
(recompression) - flow is sonic at the vena contracta. 


Pac >P = Pree (5) 


Under such a condition the fluid speed in vena 
contracta reaches the sound speed and supersonic impact 
waves arise downstream. A loud noise is given out, due to 
the fact that the sound velocity is reached and other 
complex aerodynamic disturbances are generated. 


Regime III — No recompression but with flow shear 
mechanism - flow is sonic at the vena contracta but 
supersonic after, resulting in shock waves. There is no 
isentropic recompression. 


Prec > Pe = Pop (6) 


where p2p is the valve outlet absolute pressure at break 
point. 


n pu G J 7 (7) 
and a is the recovery correction factor, given by: 

Prec (8) 
Poc 





r= 


Regime IV — Shock cell turbulent flow interaction - 
flow is sonic at the vena contracta but supersonic after. 
Shock interaction dominates noise. There is no isentropic 
recompression. 


Ps > P = Pace (9) 


where pzce is the valve outlet absolute pressure where 
region of constant acoustical efficiency begins. 


Regime V — Shock cell turbulent flow interaction - 
begins when downstream pressure drops to pzce and is 
where acoustical efficiency becomes constant. 


Pce > Pe (10) 


The velocity at the vena contracta is sonic. Further 
reductions in outlet pressure will not increase the noise 
level. 
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3.3 Jet velocity 


(11) 





where pı is the upstream specific mass of gas. 


3.4 Stream power 


The noise of interest is generated by the valve in and 
downstream of the vena contracta. If the total power 
dissipated by throttling at the vena contracta can be 
calculated, then the fraction that is noise power can be 
determined. The noise from a freely expanding jet of air 
would be equal to some acoustical efficiency factor 
multiplied by the power of the jet Wm. Equation (12) is 
valid for Regime I while equation (13) is valid for Regime 
Il, MI, IV and V. 


N {tien 
Win = Ton Ove (watts), (12) 


FẸ 
-_ 


Ww, = Teel (watts), (13) 


where mn is the mass flow and Cycc is speed of sound in the 
vena contracta at critical flow conditions, given by: 


(14) 





where R is the universal gas constant (8314 J/kmol.K), M 
is the molecular mass of gas and Tvcc is the vena contracta 
absolute temperature at critical flow conditions. 


2% (15) 


p E = 
ME pir 


3.5 Acoustical efficiency 


The acoustic efficiency factor is not constant; it is 
approximately equal to 1.10% at Mach 1 (Singleton, 1999). 
It considers the flow conditions applied through the valve 
to determine the noise generation mechanism in the valve. 
For the five defined regimes, there is dependence on the 
relationship of vena pressure contracts and downstream 
pressure. For each of these regimes, an acoustic efficiency 
is defined and calculated. This acoustic efficiency 
establishes the fraction of the total energy flow. When 
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designing a regulating valve, lower acoustic efficiency is 
desirable. For the 5 different regimes, consider: 


y= (10-*)M,,."* (16) 
Hs — a — (o-ta Ere (17) 
ny = ns = 007) (4E) (18) 


where Myc is the Mach number at the vena contracta, and 
Mj is the freely expanded jet Mach number. 


Mye = 2 (19) 
"Eve 
(20) 
M; = 





3.6 Sound power 


The sound power of a source is the total acoustic 
energy radiated by the source per unit of time. Given in 





watts. 
Wor = Niy Wn RK” (21) 
Waz = Nair Mins e] Ga 
Was = aniy Ving (23) 
Wis = Nat Wine (24) 
Was = Msi Wns (25) 


3.7 Peak frequency 


Because the transmission loss through the pipe walls is 
affected by the frequency of the sound waves it is 
experiencing it necessary to calculate the peak frequency 
of the noise generated by the valve. The peak frequency is 
given in Hertz. 


p Se (26) 
fı Di 
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| op _ Mj (27) 
D.5eyce 
fou = fox = —2 (28) 


Coo oe 


3.8 Internal sound pressure 


Sound power and sound pressure are customarily 
expressed using the logarithmic scale known as the decibel 
scale. This scale relates the quantity logarithmically to 
some standard reference. Thus, the internal sound pressure 
level, Lpi, referenced to 2 x 10° Pa for sound pressure is 
calculated in dB from the following: 


_ any. [E2107 WE p20: 2 
Lpi = 10logi n] = 


where p2 is the gas density and c2 is the speed of sound, 
both at downstream conditions. Dj; is the internal 


downstream pipe diameter. 


m =p |E] (30) 


n- PR GD 
= y M 


3.9 Pipe transmission loss 


The pipe vibration mode, for the objective of this 
prediction method is determined from a simplified model, 
considering the noise source peak frequency and the 
Natural tube 
frequencies are functions of tube diameter, wall thickness 


natural frequencies from the tube. 


and density. Transmission loss according to the model 
used by the IEC standard is based on the work of 
Fagerlund and Chow (Singleton, 1999). 


Ring frequency (f;) - as its name implies, this is the 
frequency at which a ring section of the pipe vibrates 
naturally circumferentially. Frequency that has a 
wavelength exactly equal to the circumference of the tube, 
favoring the production of a resonant stress around this 
circumference. Mathematically: 

a (32) 


a ED; 
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External coincidence frequency (f) - the frequency at 
which the external acoustic wave speed is equal to the 
velocity of a flexural wave is the pipe wall. To the speeds 
of sound of 5000 m/s and 343 m/s, in the steel and in the 
air, respectively (t = pipe wall thickness): 


f — 12.974 22 
fa : (33) 





First internal coincidence frequency (fo) - for 
frequencies higher than the cutoff frequency (wavelengths 
shorter) the acoustic pressure waves are able to travel in 
transverse direction rebounding off the pipe walls as they 
travel down the pipe. This frequency has the lowest 
natural frequency of the pipe wall. It produces a 
longitudinal flexural wave that spirals along the length of 


the pipe. 


O á 


Cutoff frequency (f.e) - frequency that is significant 
because at the cutoff frequency and below, the 
wavelengths are too long to reflect off the internal pipe 
wall, making them incapable of vibrating the pipe. 

f= ETA (35) 


Di 


3.10 Transmission loss across the pipe wall 


The only way in which noise within the downstream 
pipe can be experienced outside the pipe is through 
vibration of the pipe walls. If there is no vibration, no 
sound will be transmitted. Transmission loss through the 
pipe wall can be calculated by the equation: 


TL = 10log,, | 7.61077) (=) majl 
buy |= æ= 
x oi Sa Y} 


tf 





where Gx and Gy are the frequency factors (according to 
“Table 2”), pa is the atmospheric pressure outside pipe and 
ps is the standard atmospheric pressure. 


Table.2: Frequency factors (Lipták, 2006 - adapted) 
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fp < fo fp > fo 


= l G.=1>f,>f 
6 =F 2h <I Ip =f 


Ip | 
Gy = > h< fa 


Gy sia Ío = iy 


Gy =1=> fo = fg 


3.11 Downstream pipe velocity correction factor 


Downstream pipe velocity correction factor can be 
calculated by the equation: 


Ly = 16logy9 (——) (37) 


1—Ms 


where M2 is the Mach number in downstream pipe. 


M, = =m (38) 


-wD Pata 


3.12 Weighted sound pressure level 1 m from pipe wall 


The recognized observation point for noise predictions 
is 1 m downstream from the valve and Im from the pipe 
wall (the pipeline is considered as a line source). The 
sound pressure level is adjusted, with a 5 dB correction to 
account for all frequency peaks, for this location by: 


; | | fD ++ 
Lpaeam = 5+ Lpi +TL + Ly — log (EE) 69 


IV. COMPUTATIONAL ROUTINE 


The distributed natural gas comes from different 
production points, configuring variations in the average 
chemical composition of the gas, a condition that affects, 
even in small proportions, the calculations performed in 
the context of the present work. Thus, when proposing the 
corresponding computational routine, it was decided to 
allow the end user to adjust the partial percentages of the 
various independent components that might be part of the 
final chemical composition of the gas considered. As an 
example, the average chemical composition of the gas 
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used in the reference calculations of the present work is 
reproduced in “Table 3”. These values are compatible with 
the average chemical composition of natural gas from 
Bolivia. 


The computational routine developed has the main 
objective of facilitating the effective evaluations and the 
necessary calculations, by the operations engineer 
involved. From simplified input data it is possible to obtain 
answers about the predominant flow rate type as well as 
the calculation of the noise generated by the pressure 
regulating valve. 


Table.3: Average chemical composition of the natural 
gas considered 





Component Mole fraction Molar mass 

% kg/kmol 

CH4 16.043 
C2H6 30.069 
C3Hs 44.096 
n-C4Hio0 58.123 
1-C4H 0 58.123 
n-CsH)2 72.151 
1-CsH)2 72.151 
n-Ce6Hi4 86.178 
N2 28.013 
CO? 44.010 
Natural gas 18.374 


In sequence, as an example, the result of a computational 
routine performed. 


***Program for Regulators Valves Noise Calculation - 
according International Standard IEC 60534-8-3*** 


***Preliminary data entry ****## eee RRR RRR EE 
Chemical composition of natural gas analyzed: 
CHa: 89.01 

C2H6: 5.93 

C3Hs: 1.85 

nC4Hı0: 0.42 

1C4H jo: 0.31 

nCsHj2: 0.11 

iCsHy2: 0.08 

nC6H14: 0.08 

O2: 0 
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No: 0.67 

CO: 1.54 

total percentage: 1.00 

critical pressure of natural gas - Pc (MPa): 4.6316 
critical temperature of natural gas - Tc (K): 204.4777 


specific heat the constant pressure of natural gas - Cp: 
2177.8405 


specific heat at constant volume of natural gas - Cv: 
1690.4426 


ratio between specific heats of the natural gas - (Cp/Cv): 
1.2883 


molecular mass of the natural gas - M (kg/kmol): 18.3745 
enter the value of Fa:0.31 

enter the flow coefficient of the valve - Cy:152 

specific mass of the fluid upstream - pı (kg/m°):10.3642 
jet diameter - D; (m): 0.016679 


*** Calculation of pressure and pressure 
TALES * >% 5K EE k k A k k k ote k a k oS k k sO ie k k k K k eho ai k a ris K a ie aie a K K 


pressure of "vena contract - Pvc (MPa): 1.5123 


pressure of "vena contract" at critical flow conditions - 
Pvcc (MPa): 1.3696 


downstream pressure at the critical pressure drop (sonic 
flow) - Pac (MPa): 1.5844 


break point pressure (shocks turbulent) - Pog (MPa): 
0.93234 


factor correction: 0.86444 


downstream pressure (where acoustical efficiency 
becomes constant) - P2ce (MPa): 0.13146 


Subsonic flow; isentropic recompression; turbulent shear 
noise - Regime I 


gas velocity for the "vena contracta" region - Uvc (m/s): 
478.8963 


stream power for the "vena contracta" region - Wm 
(Watts): 1043504.4265 


absolute temperature for the "vena contracta" region - Tcv 
(K): 476.3826 


speed of sound for the "vena contracta" region - cyc (m/s): 
526.9724 


Mach number for the "vena contracta" region - Myc: 
0.90877 


acoustical efficiency factor - nı: 7.0865.10 
sound power - W, (watts): 14.9745 
peak frequency - fp (Hz): 5742.604 
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*-*E-ENT OSE calculations ** % % 7 2 7 2 2 2 2 2 ie HI ae 2 ae K k k K k k kk K K 


downstream mass density - p2 (kg/m°): 7.0477 
downstream sonic velocity - c2 (m/s): 550.6218 

Mach number at valve outlet - Mo: 0.28924 

internal sound pressure level - Lpi (dB): 162.522 

ring frequency - fr (Hz): 15603.4258 

external coincidence frequency - fg (Hz): 2154.9834 
first internal coincidence frequency - fo (Hz): 6262.0888 
cutoff frequency - fe (Hz): 3163.376 

frequency factor - Gx: 0.38479 

frequency factor - Gy: 1 

transmission loss across the pipe wall - TL: -51.4464 
downstream pipe velocity correction factor - Lg: 2.3504 


weighted sound pressure level Im from pipe wall — 


V. CONCLUSION 


This paper presents a simplified and fast proposal for 
the determination of aerodynamic noise in pressure 
regulating valves used in piped natural gas distribution 
networks. Such a proposition, focused on the noise 
prediction standard IEC 60534-8-4, accomplished the 
difficult task of producing reasonably accurate answers to 
a more complex problem entirely by theoretical means. 
The possibility of variation in the input parameters of the 
problem also allows the professional involved, a more 
accurate comparative analysis from variable operational 
data. 


It should also be noted that the IEC 60534-8-4 - 
Prediction of noise generated by hydrodynamic flow - is 
the most reliable method available at the moment. 
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